Realization of logically labeled effective pure states for bulk quantum computation 
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We report the first use of "logical labeling" to perform a 
quantum computation with a room-temperature bulk system. 
This method entails the selection of a subsystem which be- 
haves as if it were at zero temperature — except for a de- 
crease in signal strength — conditioned upon the state of 
the remaining system. No averaging over differently prepared 
molecules is required. In order to test this concept, we exe- 
cute a quantum search algorithm in a subspace of two nuclear 
spins, labeled by a third spin, using solution nuclear magnetic 
resonance (NMR), and employing a novel choice of reference 
frame to uncouple nuclei. 

The extraordinary potential of quantum algo- 
rithms to solve certain problems in fewer steps than 
the best known classical algorithms, has recently been 
highlighted by the experimental demonstration of ele- 
mentary quantum computations |p|-pT|. These exper- 
iments made use of nuclear spins, subject to a mag- 
netic field, as quantum bits (qubits: 2-level quantum sys- 
tems). Even though traditional quantum algorithms as- 
sume that all qubits are necessarily initialized in a pure 
state, typically the ground state, the experiments were 
carried out at room temperature and the spins were in a 
highly mixed state. 

The use of room temperature nuclear spins for quan- 
tum computing relies on the ideas of Gershenfeld and 
Chuang [Q and Cory, Fahmy and Havel jl3| to convert 
the mixed thermally equilibrated state into an "effective 
pure" state, at the expense of a signal decay, exponential 
the number of qubits. The density matrix of such a state 
is diagonal and proportional to the identity matrix, ex- 
cept that one of the diagonal elements is different from 
the others. Because states proportional to the identity 
matrix are not observable in NMR spectroscopy and do 
not transform under unitary operations, an effective pure 
state gives the signal and has the dynamical behavior of 
a pure state. Until now, the signal of an effective pure 
state has been extracted by adding the signals from dif- 
ferently prepared molecules, either in different locations 
using magnetic field gradients [ fl3[ or by time-sequential 
experiments E3. 

Logical labeling offers a fundamentally different ap- 
proach for creating effective pure states: a submanifold 
spanned by 2 k states is selected from the Hilbert space 
of n coupled spins, by conditioning upon the state of 
n — k ancilla spins. Gershenfeld and Chuang jl2| pre- 
dicted that the structure in the density matrix at equi- 
librium can be exploited to make this subsystem behave 



as if it were in a pure state, even if the entire system is 
at room temperature. This concept of embedding was 
previously used to observe Berry's phase in NMR spec- 
troscopy EEL 

In this Letter, we report the experimental implemen- 
tation of logical labeling by embedding an effective pure 
2-qubit subspace within a 3-spin molecule. Control over 
the dynamical behavior of the two logically labeled qubits 
is facilitated by the "uncoupling frame" , a novel method 
to remove the effect of undesired couplings on the spin 
dynamics within a subspace. We execute Grover's quan- 
tum search algorithm as a test of logical labeling and 
the uncoupling frame, and study the preservation of the 
effective pure state after many operations. 

The Hamiltonian of weakly coupled nuclear spins in 
a molecule in solution and subject to a strong magnetic 
field, is well approximated by |1(|] (K = 1) 

H = — ^2 u ilzi + ^ 2irJijI zl I z j + Henv, (1) 

i i<j 

where 2I Z i is the z Pauli operator and uii is the Larmour 
frequency of spin i. Jij is the strength of the scalar spin- 
spin coupling and Tienv represents coupling to the envi- 
ronment, which causes decoherence. For a 3-spin system 
at typical magnetic field strengths, the Boltzmann distri- 
bution describing the thermal equilibrium populations of 
the states {|000), |001), |0I0), |011), |100), 1 101), |110), 
| 111)} (where and 1 represent a spin in the ground 
and excited state respectively) is dominated by the Zee- 
man energy shift ±7iWi/2 and, for homonuclear systems 
(u>i «w,t = l,2,3), is well approximated by 

^■{1+ 3a, 1+a, 1+a, 1 — a, 1+a, 1— a, 1— a, 1— 3a}, (2) 

where a = Tiuj /2kgT << 1, with kgT the thermal energy. 
The thermal equilibrium state is thus highly mixed. How- 
ever, the observable signal from the subspace spanned by 
the states |000>, |011), 1 101) and 1 110) is the same as the 
signal from a system in the state |000) , to within a scaling 
factor. This subspace is thus in the effective pure state 
1 000). 

In order to simplify subsequent logical operations and 
to separate the signals of the effective pure subspace and 
its complement, the populations can be rearranged by a 
sequence of 1 and 2-qubit unitary operations to obtain 

-_- { 1+ 3a, 1— a, 1— a, 1— a, 1+ a, 1+ a, 1+a, 1— 3a }. (3) 
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Now the subspace {|000), |001), |010), |011)} is in an 
effective pure state. This subspace corresponds to the 
spins B and C conditioned on or labeled by the state 
of the first spin, A, being |0) (we will call this the | 0) A 
subspace). This procedure is called logical labeling |l2|] 
and, combined with removal of coupling to spin A for the 
remainder of the pulse sequence, allows 2-qubit quantum 
computations involving only the spins B and C . 

We selected bromotrifluoroethylene dissolved in 
deuterated acetone (10 mol%) as the central molecule 
in our experiments, because the spin-1/2 19 F nuclei have 
large J-couplings and chemical shifts, as well as long co- 
herence times, which makes it suitable for quantum com- 
putation p|Jl^|. The 12 C nuclei are non-magnetic and 
the interaction of the spin-3/2 Br nucleus with the fluo- 
rine spins is averaged out due to fast Br relaxation jl6) . 
The 19 F Larmour frequencies are » 470 MHz (at 11.7 T) 
and the spectrum is first order, consisting of three well- 
separated quadruplets (as in Fig. ||), with uja — ujb/^tt ~ 
13.2 kHz and luq — lua/2tt « 9.5 kHz. The coupling con- 
stants are measured to be Jab = —122.1 Hz, J ac = 75.0 
Hz and Jbc — 53.8 Hz (see also @). Single-spin ro- 
tations are realized by applying a radio-frequency (RF) 
magnetic field, oscillating in the x — y plane. To simul- 
taneously address the four lines in one quadruplet with- 
out affecting the other two quadruplets, the envelope of 
the RF pulses is Gaussian-shaped (more complex shapes 
could improve the results) and the RF power is adjusted 
to obtain pulses of ss 300/is. Experiments were performed 
at IBM using an Oxford Instruments wide-bore magnet 
and a Varian UNITY Inova spectrometer with a Nalorac 
triple resonance (H-F-X) probe. 

The logical labeling step is implemented as follows. 
Comparison of Eq. ^ with Eq. |^ shows that the popu- 
lations of the states |001) <-> |101) and |010) «-> |110) 
must be interchanged, while the remaining four popula- 
tions must be unaffected. This requires a CNOT^ and 
a CNOTca, where CNOT„ represents a controlled-NOT 
which flips j if and only if i is in |1). Furthermore, the 
two CNOT gates commute and can thus be executed si- 
multaneously, as shown in Fig. [l]. 
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FIG. 1. The logical labeling pulse sequence. Y and Y de- 
note a 90° rotation about the y and — y axis (right hand rule) . 
X 2 are 180° rotations. 

This sequence implements the CNOT gates only up 
to single-spin Z rotations, which suffices for a diagonal 
initial state as used here. Furthermore, the I z b, IzC an d 
2I z bI z c terms in the Hamiltonian have no effect, since 
B and C remain along ±z. I z a can be ignored because 
the pulses were applied in a reference frame in resonance 



with each spin. Fig. || shows the measured populations 
of the eight basis states, before and after the sequence of 
Fig. |l|. The results agree with the theoretical predictions 
of Eqs. | and |. 
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000 001 010 011 100 101 110 111 000 001 010 011 100 101 110 111 
FIG. 2. Experimentally determined populations (in ar- 
bitrary units, and relative to the average) of the states 
|000), . . . , 1 1 1 1 } (Left) in thermal equilibrium and (Right) af- 
ter logical labeling. The populations were determined by par- 
tial state tomography [ po[ . 

The effective pure state prepared with logical labeling 
must remain pure throughout any subsequent computa- 
tion. This requires that while a computation is carried 
out using B and C, A must "do nothing" , which is non- 
trivial in a system of coupled spins ]is| : the effect of 
Jab and Jac must be removed. This could be done by 
using two refocusing X\ pulses during every logical op- 
eration between B and C. We have devised a different 
method, which exploits the fact that it suffices to remove 
the effect of Jab and Jac within the |0) A subspace. This 
uncoupling frame method requires no pulses at all and is 
described in Fig. ||. 
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FIG. 3. Spectrum of the transitions of B, with the states of 
the other two spins as indicated. With respect to a reference 
frame rotating at lob/^tt —Jab/2 (indicated by an arrow), 
spins B which see a spin A in |0), evolve under Jbc only and 
are thus uncoupled from A. Jab does affect _B's evolution in 
the |1) A subspace, but the signal of this subspace does not in- 
terfere with that of the |0) A subspace. Similarly, C's rotating 
frame must be moved to ujc/^tt — Jac /2. A separate channel 
was used for B and C. We note that the uncoupling frame 
is somewhat related to selective decoupling a technique 
to determine the relative sign of J-couplings by moving the 
reference frame of one nucleus (not several, as here) to the 
center of a submanifold. 

Mathematically, the transformation of spin i's state 
from the laboratory frame to a reference frame rotating 
at oji + Aj (with the resonance off-set frequency), is 
described by \ip) = U\ip') — exp[z(cji + Ai)I zi t]\ip'} , where 
IV^') is the state in the new reference frame. After sub- 
stituting by U\ip'} in the Schroedinger equation, and 
using the commutation rules between I x , I y and I z , it 
follows that H' = H + h(uJi + Ai)I zi Then, choosing 
A B = -ttJab and A c = -kJac gives 



H' = 2n[(l/2 + I zA ) JbcLbI 



zC 
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+ (1/2 — I z a) {JbcIzbIzC - JabIzB — JacIzc) ] , 

where 1 is the identity matrix. The first (second) term in 
the expression of TC acts exclusively on the |0) A (ll)^) 
subspace. Any state within the |0) A subspace evolves 
only under JbcIzbIzC and will remain within the 10)^ 
subspace. Within the |0)^ subspace and using the uncou- 
pling reference frame, B and C are thus uncoupled from 
A. We experimentally confirm this by reconstructing the 
deviation density matrix p"2"| of the 3-spin system after 
creating the state ( 1 00) + |ll))/-\/2 in the |0) A subspace. 
Fig. U shows the |0) A subsystem in the predicted effective 
pure state and uncoupled from the subspace. 




FIG. 4. Normalized experimental deviation density matrix 
(with the diagonal shifted to obtain unit trace for the effective 
pure state), shown in absolute value. The entries in the second 
quadrant are very small, which means that the |0} A and |1} 4 
subspaces are uncoupled. The deviation density matrix is 
measured using quantum state tomography H ; a series of 27 
consecutive experiments with different sets of read-out pulses 
(see below) permits to determine all the components of the 
deviation from identity of the density matrix. The relative 
error in the state is || p exp — Pth ||/|| Pth || = 19 %. 

The motivation for creating effective pure states is to 
provide a suitable initial state for quantum computation, 
and as a test of this, we implement Grover's quantum 
search algorithm || on the two logically labeled spins, 
using the uncoupling frame. With this algorithm, one can 
find the unique but unknown xq which satisfies f(xo) = 1, 
in 0{s/~N) queries of f(x), where N is the number of 
possible entries x. Searching on a classical machine would 
require O(N) queries of f(x). For TV = 4, implemented 
here, a quantum computer finds xq in just one query, 
compared to on average 2.25 queries classically. 

We use a pulse sequence similar to the one described 
in ||. First, YbYq rotates both spins from 1 00) into 
(| 00) + |01) + 1 10) + |ll))/2, an equal superposition of the 
four possible inputs. The amplitude of the \xq) term is 
then amplified in two steps ||. First, one of four func- 
tions f Xo ( x ) is evaluated, flipping the sign of the \xq) 



term. This is done by one of four conditional phase flips 
Y B Y C - $ B 6 C - Y B Yc - 1/2J BC , where $ = X for / 00 
and /io and $ = X for f 01 and f n . 9 = X for / 00 
and foi and = X for /io and fn. Second, inversion 
about the average || is implemented by a Hadamard 
gate (a tt rotation about an axis halfway between x and 
z) on both spins, followed by the conditional phase flip 
corresponding to /oo, and another Hadamard gate. The 
pulse sequence for this inversion step can be reduced to 
XbXc — YbYq — \/2Jbc ~ YbYq- The entire sequence 
for Grover's algorithm takes approximately 20 ms and the 
labeling step takes about 7 ms. The coherence time for 
the three 19 F spins, expressed as the measured transverse 
relaxation time constant T2 ~ 4-8 s, is sufficiently long 
for coherence to be maintained throughout the labeling 
and computation operations. 

The theoretical prediction is that the state of BC after 
completion of the algorithm is the effective pure state 
\x(s), which can be determined by a measurement of the 
spin states of B and C. This is done using a read-out 
pulse for each spin, which rotates the spin back into the 
x — y plane so that it generates a rotating magnetic field 
and induces an oscillating voltage in a pick-up coil. The 
voltage is recorded with a phase sensitive detector, and 
Fourier transformed. The phase of the receiver is set 
to obtain a positive (negative) absorption line for a spin 
that was |0) (|1)) before the read-out pulse, and is used 
to determine |a;o) (see Fig. ||). The experimental spectra 
(Fig. H Left) as well as the deviation density matrices of 
the logically labeled subspace before, during and after the 
computation, confirm that the state remains an effective 
pure state throughout the computation and that the final 
state is |cco). 

An interesting question is how many logical operations 
can be executed while preserving the effective pure char- 
acter of the spins, and further, how quickly errors ac- 
cumulate during longer pulse sequences. We study this 
by iterating the conditional flip and inversion steps in 
Grover's algorithm, which ideally gives rise to a periodic 
pattern: for N = 4, the amplitude of the xq term is ex- 
pected to be 1 after I iteration, and again after 4, 7, . . . 
iterations M. Fig. || demonstrates the expected periodic 
behavior in the output state in experiments with up to 
37 iterations, which requires 448 pulses and takes about 
700 ms. In addition to decoherence, errors mainly arise 
from imperfections in the pulses and coupled evolution 
during the 300 /is pulses (reduction of t\ and T2 partially 
compensated for this effect). Spectra of the quality of 
Fig. |^ were obtained by choosing a particular implemen- 
tation of the composite z rotation in the phase flip step, 
such that the errors it introduces partially cancel with the 
errors of the inversion step. For example, while YXY, 
YX Y, XYX and XYX are all mathematically equiva- 
lent, the errors may in practice add up or cancel out with 
the errors from previous or subsequent pulses. Clearly, 
a general optimization procedure will be very helpful for 
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designing effective pulse sequences in future experiments 
involving more qubits. 
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FIG. 5. Real part of experimental spectra (frequencies in 
the uncoupling reference frame) for spin B (Top) and C (Bot- 
tom), after executing Grover's algorithm 1 (Left), 19 (Center) 
and 37 (Right) times, with |a;o) = 1 00) . The |0) A subspace 
corresponds to the spectral lines at ±Jsc/2 = ±26.9 Hz. Ide- 
ally, the line at -26.9 Hz is positive and absorptive, with unit 
amplitude, while the line at +26.9 Hz is zero. Even after 37 
iterations, xq can be unambiguously determined. 

We attribute the preservation of effective pure states 
during long pulse sequences in part to the use of the 
uncoupling frame, which provides an elegant and sim- 
ple alternative to refocusing schemes involving 7r-pulses. 
Multiple couplings with ancillae spins can be neutralized 
simply by moving the carrier frequencies of the compu- 
tation spins by the appropriate ^±J/2. In contrast, 
refocusing pulses, would have to be applied during every 
single evolution interval and their complexity rapidly in- 
creases as more J-couplings are to be refocused Jl8| . This 
technique may find application in future experiments us- 
ing logical labeling, as well as in quantum error detection 
experiments |2l| , where the computation must only pro- 
ceed within the subspace labeled error-free by the ancil- 
lae. However, refocusing pulses are still required when a 
coupling must be removed over an entire system rather 
than in a subspace only. 

The present experiments demonstrate that a fc-qubit 
room temperature system can behave as if it were very 
cold, except that the signal strength decreases exponen- 
tially with the number of qubits, when it is properly em- 
bedded in an n-spin system. The subspace dimension is 
limited by the number of equally populated states, which 
is Cn^ 2 — nl/(n/2)l 2 in a homonuclear system, giving 
k = log 2 (l + Cl l J 2 ). Thus for large n, k/n tends to 1 
(n=40 for fc=37). For heteronuclear spin systems, the 
analysis is more complex and k/n is generally smaller, 
but the number k of cold qubits that logical labeling can 
extract from n hot spins, as well as the number of oper- 
ations required for logical labeling, still scale favorably. 
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